Detecting dark matter in electromagnetic field penetration experiments 
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Dark matter in the form of particles from a hidden mirror sector has recently been proposed as 
an explanation for the DAMA annual modulation signal. Here one assumes that there exists a small 
mixing between photons and mirror photons. We show that dark matter with this property can 
also be detected in electromagnetic field penetration experiments. Such experiments can be used to 
measure the speed and direction of the dark matter halo wind, the local density, the temperature, 
and the strength of the photon-mirror photon mixing interaction. An additional result would be a 
significant improvement of the upper limit on the photon mass. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

Dark matter may couple to the electromagnetic field. 
Although charged dark matter is strongly constrained , 
particles with magnetic and/or electric dipole moments 
Q, axions Q and millicharged dark matter Q are vi- 
able dark matter candidates. Another dark matter candi- 
date that can couple to the electromagnetic field is mirror 
matter. Mirror matter is predicted by extensions of the 
Standard Model that restore parity symmetry 0, @ ■ In 
these models the Standard Model is appended by a hid- 
den sector which consists of a parity transformed copy of 
the Standard Model. Parity can either be spontaneously 
broken or remain unbroken depending on the Higgs po- 
tential ■ In this article we shall consider the unbroken 
case where each Standard Model particle has a mirror 
partner of exactly the same mass. The doubling of the 
number of relativistic degrees of freedom during Big Bang 
Nucleosynthesis (BBN) requires that the temperature in 
the mirror sector be lower than in the ordinary sector by 
at least a factor of 2. It has been shown that this can be 
a typical situation after the inflationary epoch 0]. Vari- 
ous aspects of mirror dark matter are discussed in recent 
reviews U E3 and in Refs. H, EH. See Ref. [H for 
a historical review. Non-minimal mirror matter models 
consisting of multiple copies of the Standard Model are 
discussed in Ref. |13| . 

Mirror matter would be expected to form compact ob- 
jects. Such objects could be responsible for the large 
number of observed MACHOs with masses M ~ O.5M 
[l4| . close-in extrasolar planets 0] and isolated plan- 
etary mass objects 0|. Non-gravitational interactions 
between mirror particles and ordinary particles could 
provide for additional signatures of the mirror world. 
Such interactions are strongly constrained by renormal- 
izability, gauge invariance and parity symmetry. In Refs. 
[6lll7| all such interactions were identified. In this paper 
we will focus on one of the possible interactions that sat- 
isfies these constraints, the photon-mirror photon kinetic 



mixing 



interaction H E E E 13 : 



C = ~ 2 F-F>, 



(1) 



where F (F 1 ) is the field strength tensor for electromag- 
netism (mirror electromagnetism). This term causes a 
mirror electric ch arge o f q appear to be an ordinary elec- 
tric charge of eq 0, [T^, • The photon- mirror photon 
kinetic mixing interaction has been invoked to explain the 
DAMA annual modulation signal [2l| , gamma ray bursts 
|2^ . problems with the central engine of supernova ex- 
plosions and the 511 KeV gamma-rays from the galac- 
tic bulge [22^ . the cosmic coincidence problem |23j . the 
pioneer spacecraft anomaly 0> Epl- the low number of 
small craters on the asteroid Eros [25j and anomalous im- 
pact events |25l l2ti| . The possible role of neutrino- mirror 
neutrino oscillations in gamma ray bursts is discussed 
in Refs. |27j and neutron-mirror neutron oscillations has 
been postulated as an explanation for ultra high energy 
cosmic rays [28). 

The photon-mirror photon kinetic mixing interaction 
JTJ can affect BBN. The agreement of standard BBN the- 
ory with observations implies a limit ofe < 3x 10~ 8 [29| . 
Another important effect of the kinetic mixing term is 
that it induces orthopositronium-mirror orthopositron- 
ium oscillations |19|| which can be detected in vacuum 
experiments [3(J. The oscillations lead to a higher ef- 
fective decay rate compared to the QED prediction, be- 
cause the mirror photons from the decays of mirror or- 
thopositronium are not detected. Because collisions be- 
tween orthopositronium and gas molecules will inhibit 
the oscillations, it is possible to distinguish this effect 
from possible invisible decay modes of orthopositronium. 
In a vacuum cavity experiment performed in 1990 |3l| 
a faster decay rate than the QED prediction was seen. 
This was also in conflict with the results other experi- 
ments in which the cavity was not vacuum and which 
were consistent with the QED prediction. A value for |e| 
of about 1CP 6 is needed to explain this anomal y 1321 . In 
2003 a new vacuum experiment was carried out [33j , and 
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this time the anomaly was not seen. This latest exper- 
iment implies that e<5xl0~ 7 [j|. A more sensitive 
experiment is currently being planned |30j |. 

If |e| > 10~ 10 then macroscopic mirror fragments can 
remain close to the Earth's surface and be detected using 
centrifuges Q. It would also be possible to locate such 
fragments using the fact that these fragments draw in 
heat from their surroundings, radiating it away as mir- 
ror photons and thereby causing an observable cooling 
effect |35| 1 . One would also expect observable effects 
from meteorites [2^, and micrometeorites [^fj. If 
l e l ^ then heavy atoms can bind to heavy mir- 

ror atoms [33]. Such bound states could be detected as 
anomalously heavy isotopes. Other possible effects of 
bound states were considered in Ref. |3Sj . 

Another detectable feature of the mirror world is pro- 
vided bv mirror ions in the galactic halo. As explained in 
Refs. |9ll2l|. most of the mirror dark matter in the galac- 
tic halo should be in the form of a plasma. It is expected 
that this plasma consists mainly of mirror helium and 
heavier elements, because nucleosynthesis in the mirror 
sector is expected to yield a large abundance of mirror 
helium Assuming that the dark matter halo of 

our galaxy consists of a mirror helium plasma, hydro- 
static equilibrium then implies that the temperature of 
this plasma is about 330 eV @. Such a hot plasma will 
cool itself down rapidly and collapse into a disk, unless a 
heat source of about 10 36 — 10 37 Watt exists to keep the 
mirror plasma hot. Supernovae of mirror stars at a rate 
of about 1 per year could provide for such a heat source 
. Such a high rate of mirror supernovae could be the 
result of the exp ected high mirror helium abundance of 
mirror stars [39j. 

The galactic mirror matter plasma may already have 
been detected in the DAMA/Nal experiment 0, liof . 
The DAMA/Nal experiment is one of about 20 experi- 
ments in which one tries to detect recoils of nuclei caused 
by collisions with dark matter particles. The positive re- 
sults of the DAMA/Nal experiment are difficult to recon- 
cile with the negative results of other experiments, such 
as e.g. CDMS II ^l| if one assumes that dark matter 
consists of neutralinos. The crucial difference between 
neutralinos and mirror ions is that neutralinos are ex- 
pected to be much heavier. Unlike other detectors, the 
DAMA detector is capable of detecting recoils caused by 
mirror ions because it contains the light element sodium. 
The mirror world explanation for the DAMA results as- 
sumes that the detected recoils are caused by mirror oxy- 
gen ions or heavier ions colliding with sodium ions in the 
DAMA detector. Assuming that mirror oxygen is re- 
sponsible for the DAMA signal, the value for |e| needs to 
be about 5 x 10~ 9 -\/0.1/£, where £ is the relative abun- 



dance by mass of mirror oxygen |21j . Other explanations 
for the DAMA results include light WIMPs Hg, WIMPs 
with spin dependent interactions 0] , inelastic dark mat- 
ter [44|, four th g eneration neutrinos [43, self-interacting 
dark matter [46j , and pseudoscalar and scalar dark mat- 
ter H3. 

In this article we will explore another way of detecting 
the mirror plasma. The photon-mirror photon kinetic 
mixing term couples the mirror plasma to ordinary 
electromagnetic fields. This means that in very strong 
electric fields the plasma will become slightly polarized. 
Mirror charges will be induced which one can attempt to 
detect inside superconducting cavities. Strong magnetic 
fields should also give rise to polarization effects, because 
the plasma is expected to have a non-zero velocity, Vh, 
w.r.t. the Earth. By comparing the results of experiments 
using magnetic and electric fields one can thus measure 
both the magnitude and the direction of The velocity 
of the plasma also makes the polarization effect in purely 
electric fields slightly anisotropic. We will show that this 
effect can be used to determine |e|, the temperature and 
the density of the mirror plasma. 



II. POLARIZING THE MIRROR PLASMA 

Unlike direct detection experiments which arc sensitive 
to the mirror ions, polarization effects are mainly caused 
by the response of the mirror electrons to the electromag- 
netic fields. The most important quantities determining 
this effect are the mirror electron density and the tem- 
perature. Assuming that the galactic mirror plasma is 
spherical, isothermal and in hydrostatic equilibrium, the 
temperature is: 



T 



1 



(2) 



where rh is the average mass of the mirror ions and the 
mirror electrons and \v m t\ ~ 220 km/s is the rotational 
speed @ . For a mirror plasma consisting of purely mirror 
helium, T w 330 eV. It is conven ient to define the electron 
thermal speed as v = y/2T/m e « 1.1 x 10 7 v /T/(330 eV) 
m/s. Boltzmann factors of exp[— m e u 2 /(2T)] can then 

be written as exp —(v/vq) . 

Most light elements except hydrogen that have a large 
abundance in the universe have a nuclear mass of about 
twice the proton mass. Assuming that this is also the 
case for mirror matter, we can relate the mirror electron 
number density, n e , to a good approximation to the total 
mirror dark matter density, p, and the contribution of 
mirror hydrogen to p which we denote by pu- 
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2m p n e - p H - (3) 



1 Since no mirror stars can be inside our solar system, we don't 
have to consider absorption of mirror photons by the mirror frag- 
ment. 



Here the nj denote the number densities of ions with mass 
m,j and m p is the proton mass. If pu can be ignored then 
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TABLE I: Properties of the mirror plasma assumed in this 
article. 



Composition 


Mirror helium 


Temperature 


T = 330 eV 


Mirror electron thermal speed 


v = 1.1 x 10 7 m/s 


Mirror electron density 


n e = 0.16 cm -3 


Plasma angular frequency 


uj p = 23 x 10 3 Hz 


Debye screening length 


A D = 340 m 


Halo wind speed 


|«h| = 220 km/s 



n e w 0.16/cm 3 if p — 0.3 GeV/cm 3 . Later in this article 
we'll see that n e can be determined experimentally. The 
hypothesis p « 0.3 GeV/cm 3 and pn ~ can thus be 
tested. More generally, upper and lower bounds for pn 
can be set. 

The plasma angular frequency for n e « 0.16/cm is 
uj p = \/n e e 2 1 (meo) ~ 23 x 10 3 Hz. The Debye screening 
length is Ad = 2~ 1 / 2 vq/uj p sa 340 m. The frequency 
of Coulomb collisions is ~ ui p log(A) /A where A which 
is proportional to the average number of particles in a 
Debye sphere: A = 47rA D n [4S|. Since Aw6x 10 13 , the 
plasma is nearly collisionless. In Table [i] we have listed 
the values for the plasma parameters that we'll use in 
this article. 

We can roughly estimate the magnitude of polariza- 
tion effects as follows. If we place an ordinary electric 
charge q in a mirror plasma then, ignoring any distor- 
tions due to the motion w.r.t. the plasma, the mirror 
electric potential will be: 



V(r) 



q 



exp 



A D 



(4) 



~4tt£o \r\ 

Consider a capacitor consisting of two infinite conduct- 
ing plates. If an electric potential difference is maintained 
between the two plates, then it follows from a straightfor- 
ward integration of that outside the capacitor, close 
to one of the plates, the induced mirror charge density 
is: 

P=-ee V^ = -f ^. (5) 

Here AV is the potential of the plate closest to the point 
where we evaluate p relative to the plate furthest from 
that point and we have assumed that the distance be- 
tween the two plates is much smaller than Ad- This 
mirror charge density will also exist inside a supercon- 
ducting cavity placed close to the capacitor and will give 
rise to electric fields in the interior of the cavity. In case 
of a cylindrical cavity with a cylindrical conductor at the 
center, the ordinary electric potential difference between 
the outer and inner conductor, AV^ av is: 



2 J2 



AV Ci 



e r. 



\2 
A D 



(6) 



Here ri and T\ are the radii of the outer and inner con- 
ductors, respectively. Electric field strengths between 
the capacitor plates of order 5 x 10 6 V/m can be eas- 
ily achieved if a material with a high dielectric strength 
such as e.g. polyethylene is used as a dielectric medium. 
Taking the radius of the cylindrical cavity r 2 ~ 1 m and 
the distance between the capacitor plates ~ 10 m, then 
values for AT4av of about 10 -15 V can be obtained if 
|e| ~ 5 x 10~ 9 (the value favored by the mirror world ex- 
planation of the DAMA results) . Static electric fields of 
this magnitude are difficult to detect. We can, however, 
let Ay oscillate at some angular frequency u>. This will 
induce an emf inside the cavity. In the next section we 
will see that the emf is approximately given by © for 
lu < Lu p and falls off to zero for ui — > oo. The emf has 
a large internal impedance caused by the small capaci- 
tance of the cavity. Suppose that the cylindrical cavity 
is curved into a toroidal shape with length I 3> r^- Such 
a cavity has a capacitance C of: 



C = 27T£ - 



nF. 



(7) 



To extract the maximum amount of power from this emf 
one can use a superconducting coil with high Q factor 
connected between the inner and outer conductors to 
cancel out the imaginary part of the impedance of the 
emf. The self inductance of the coil needs to be: 



1 /20KHz\ /30 m 

Lco . 1 = __ = 1.5l__l l_ 



log( — 

n 



H. 



(8) 

The quality factor of such a coil can be about Q ~ 10 6 
at frequencies in the KHz range |48| . The losses in the 
superconducting cavity walls are very low and can be 
ignored; Q can be ~ 10 10 for resonant cavities operat- 
ing at GHz frequencies |50| and unlike in case of non- 
superconducting cavities, the Q factor becomes higher at 
lower frequencies. The power that can be extracted from 
the emf is thus given by: 



P = "^j— (AV^v) 

AV 



5 x 10" 9 
uj \ I I 



1 / 340 m\ 4 /r., . i 



J m \ /ray 
D / ^m/ 



Q_ 

5xl0 7 V/ V20 KRzJ V 30 m J \l0 6 



lo g(^f) 



1.5 x 10" 29 W. 



(9) 



If a noise temperature of T no i SC ~ 10 K can be achieved 
then the polarization effect can be detected after an inte- 
gration time of 4fc b T noisc /P — 400 days if |e| — 5 x 10~ 9 . 

Results from this experiment can also be used to sig- 
nificantly lower the upper limit on the photon mass. Po- 
larization effects from the mirror plasma are similar to 
the effects of a photon mass of |e| /Arj on length scales 
much smaller than Ad ■ As explained in Ref. |5l| , if the 
photon acquires a mass due to the Higgs effect then the 
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sharp limits posed by the galactic vector potential are 
invalid. In that case the sharpest limit on the photon 
mass is given by the most precise experimental test of 
Coulomb's law to date [H^. which implies an upper limit 
on the photon mass of ~ 10~ 14 eV. This corresponds to 
an upper limit on |e| of 2 x 1CP 5 . Testing the mirror 
world explanation of the DAMA annual modulation sig- 
nal would lower the upper limit on the photon mass to 
5 x 10~ 9 /(340 m) w 3 x 10" 18 eV. 

Equation J^J will start to break down when w ~ w p . 
To find the polarization effect for general frequencies, we 
will set up a derivation based on the Vlasov equation. 
This will also allow for a calculation of the diurnal mod- 
ulation caused by the anisotropic mirror electron velocity 
distribution in the Earth's frame. 



III. VLASOV EQUATION 

The dynamics of a collisionless plasma is described by 
the Vlasov equation [4^ : 

df e 

2 +! ,.V/--(£ + t.xB).V„/ = 0. (10) 
at m c 

Here / is the mirror electron distribution function de- 
fined as the density of mirror electrons per unit phase- 
space volume. This has to be supplemented with the 
Maxwell equations for the mirror electromagnetic fields. 
We will ignore the generation of magnetic fields by mirror 
electron and mirror ion currents. These effects are sup- 
pressed by a factor of order \vh.\ / c relative to the main ef- 
fect. The effect of the displacement current is suppressed 
by an additional factor of order (system size) xw/c and 
will also be ignored. Within these approximations, the 
inhomogeneous Maxwell equations are: 



distribution function is: 



V • E 



so 



2m, 



fd 3 v 



' ■ 

£0 



(11) 



Here n lon is the local mirror helium ion density, /j on is 
the velocity distribution of the mirror helium ions nor- 
malized to the local density, p oxt and J ex t are the ordinary 
charge and current densities, respectively. They satisfy 
the equations: 



Pcxt — £oV ■ E ex t, 

Jext = V X -Bcxt- 

Mo 



(12) 



where E cxt and B oxt are the ordinary 'external' fields that 
are used to polarize the mirror plasma. 

Treating the external fields as a perturbation, the mir- 
ror electron distribution function can be written to first 
order as / = /q + fx where /o is the unperturbed distri- 
bution function and fx is the first order correction caused 
by the external fields. The unperturbed mirror electron 



fo(v) = n e (Trv$) 



exp 



(v - v h ) 



(13) 



Here is the average velocity of the plasma ions and 
electrons w.r.t. Earth. In case of a non-rotating plasma, 
\vb\ = u rot ~ 220 km/s. From the Vlasov equation l|10|) 
and the Maxwell equations ijTJ we obtain: 



^ + v ■ V/i - —{E + vxB)- V„/ = 0, 
at m? 



V ■ £ = eV ■ E cxt 
B = eB cxt . 



e 



d 3 vfi, 



(14) 



In appendix [S] we eliminate fx from (|14f> . If the exter- 
nal fields are proportional to sin(wi) then we find for the 
mirror charge density: 

p(x,t) =ep ext (x,t) + 

e J d 3 yA(y) [G(x - y, u, t) - G(x - y, -u, t)} , 

(15) 

where A(y) is the coefficient of sin(wt) in the 'effective 
charge': 

Pcxt (y, t) + e V • (v h x B cxt (y, t)) = A(y) sin(cjt) . (16) 
The Green's function G(x ) u> ) t) is defined as: 



2z(27r) 



K- 



\/2up Ad |fe| 



v \k\ 



(17) 



where Z' (C) is the derivative of the plasma dispersion 
function |53j and can be expressed as: 



Z'(C) - -2 + 2^Cexp(-C 2 )(erfi(C) 



Here we have defined erfi(C) as 



erfi(C) 



erf(iC) 



eft exp (t 2 ) 



(18) 



(19) 



We will now apply this to the case of a perturbation 
caused by a capacitor. Unlike in the case u> <C u> p , we 
can't ignore finite size effects here. In practice, the size of 
the capacitor will be much smaller than Ad ■ This means 
that plasma resonance effects are strongly suppressed. To 
simplify the calculations we'll take the capacitor plates 
to be disk shaped. From the estimate of the power of the 
signal at low frequencies 10, we can see that the length of 
the cavity must be of the order of ~ 30 meters. If we take 
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the cavity to be toroidal with its center on the symmetry 
axis of the capacitor, the mirror charge density is mea- 
sured at ~ 5 meters off the symmetry axis. If we choose 
the radius of the capacitor plates to be a few tens of me- 
ters then this will to a good approximation be the same 
as the charge density on the symmetry axis. Because we 
need to reserve some space between the capacitor and the 
cavity to allow for the necessary shielding and cooling ap- 
paratus, the mirror charge density should be evaluated 
a few meters from the capacitor. Evaluating the mirror 
charge density on the symmetry axis at a finite distance 
from the capacitor will thus allow for realistic estimates 
of what one can expect in experiments. 

Because |i>h| /vo » 1/50, we can obtain a good approx- 
imation for the mirror charge density by expanding it 
in this parameter to first order. If the potential of the 
plate closest to the cavity relative to the other plate is 
AV^ sin(u;i) then we find: 



p{d,u),t) = -2 — 



A x (d, w , t) + cos(0) A 2 (d, u, t) 



(20) 

Here d denotes the position where p is evaluated relative 
to the center of the closest capacitor plate and 9 is the 
angle between v-^ and d. The dimensionless functions Ai 
are given by: 



R 



/•oo 1 1 ( R \ 

Ai(d,uj,t) — 7T— 7 Im / du dflujA - — wsin(/3) )x 
2vrL J J \\ D J 



exp iu cos(/3) 



L+\d\ 



A 



D 



exp ( iu cos(P) 



Gi [ u, — ) exp(iu>t) — d { u, I exp(— itot) 

Wn / V CJ, 



(21) 



Here L is the distance between the plates, R is the radius 
of the plates and the Gi are given by: 

Gi («,/) = 



-2 M exp(^)-^/erfi(^ 

/ f9 \ i / 



iV^f - 2(u + u 3 ) exp(^) + v^F/erfi(^)'' 
G 2 (uJ) = -4u 2 exp(J^j x 
i^pnp - y/2fuexp(- 



(22) 



0.2 - 



0.4 




100 200 300 400 W P 

FIG. 1: The amplitudes ai and ai as a function of the fre- 
quency for the case R — 50 m, L — 10 m and \d\ = 4 m. 



The cos(#) factor multiplying A 2 in (|20|l is a periodic 
function with a period of one sidereal day which is de- 
termined by the Earth's motion relative to the galactic 
frame and the orientation of the capacitor. If the capac- 
itor plates are placed parallel to the Earth's surface and 
the cavity is placed below the bottom plate then we can 
write cos(#) as the inner product of the unit vector point- 
ing in the direction of — and the unit vector pointing 
from the Earth's center toward the detector: 



Jl 

2u- 



iy/^u 2 + V^(/ 2 - u 2 ) erfi(^y C os(0) = cos(<5) cos(a) + sin(<5) sin(a) cos[Q(t - *o)] 



^V2^/-2(^i + M 3 )exp(^) 



(23) 



We have numerically calculated the functions Ai as a 
function of to for the case R = 50 ru, L = 10 m and 
\d\ = 4 m. We represent these in Figs. ^ and |2 by their 
amplitudes dj and phase shifts 4>i relative to the external 
field: 



Ai(d, to, t) = a,i(d, to) sin[o;i + 4>i(u>)] 



(24) 



(25) 

where S w 43° is the angle of — Vh with the vector parallel 
to the Earth's axis, a is the angle between the Earth's 
axis and the position of the detector relative to the center 
of the Earth and can be expressed as 90° minus the geo- 
graphical latitude, fl — 2tt (sidereal day) -1 and to as 21 
hours local sidereal time so that f2(t — to) is the angle 
between the vector —v^ projected onto the Earth's equa- 
torial plane and the position of the detector relative to 
the Earth's center projected onto this plane. In appendix 
[5]we show how these values for S and to are obtained from 
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and 



fi =0i H cos(o) cos(ev) sm(02 

a x v 




400 "p 




FIG. 2: The phase shifts <f>i and 02 as a function of the fre- 
quency for the case R = 50 m, L — 10 m and |d| = 4 m. At 
zero frequency 0i = mod (27r) and 02 = 71" mod (2tt). 



the orbital parameters of the motion of the Earth around 
the Sun and the Sun around the Milky Way. 

If we insert B24JI and (|25|) in (1201) and expand to first 
order in |t>h| /vq we obtain: 



p{d, tJ,t) = -J — {h sm(uit + ipi) + 

6 2 [sin((w - O) t + <p 2 ) + sin((w + Q) t + <p a )]} , 

(26) 



where 



f2 =92 
<f3 =02 



nt . 



(29) 

(30) 
(31) 



Comparing (12(i[) with (JHJ, we see that the power of the 
signal at frequency u> is 4& 2 P w 4a^P where P is given 
by IHJ 2 . The signal at the frequencies oj±CI is about 10 4 
times weaker. 

From Ea. llfil we see that strong magnetic fields of order 
\B\ ~ 10 T will have similar effects as electric fields of or- 
der \vh \ \B\ ~ 2 x 10 6 V/m. This is similar in magnitude 
as in the case of the capacitor we have considered; experi- 
ments using large coils can thus yield similar results as ex- 
periments using large capacitors. In the case of external 
magnetic fields the leading polarization effect will depend 
on the time averaged magnitude of D[, x 5. By compar- 
ing measurements using electric and magnetic fields one 
can thus measure both the direction and the magnitude 
ofvh- The order |^h| /t>o correction of the main polariza- 
tion effect in external electric fields can be measured by 
performing experiments at different geographic latitudes. 
This thus allows one to determine vq ■ Once this parame- 
ter is fixed, the frequency profile of the polarization effect 
depends, up to an overall coefficient proportional to e 2 , 
only on uj p because Ad = 2~ 1 / 2 uo/w p . A fit to the fre- 
quency profile is thus sufficient to determine u> p , Ad and 

If Xd is much larger than the size of the capacitor or 
the coil, then it will be difficult to determine uj p and 
Ad- In that case the polarization effect becomes almost 
invariant under the rescaling: Ad — * S Ad, ~* u p /S, 
and e — * Sc. If oj p can be measured then that fixes the 
mirror electron density and, as explained in Section [HJ 
this allows one to infer the mirror dark matter density 
(see Eq. El . Effects due to mirror matter can be easily 
distinguished from effects due to a finite photon mass 
because in the latter case the frequency profile is flat and 
the phase shift is zero. 



IV. DISCUSSION 

We have shown that in strong electromagnetic fields 
mirror dark matter can become polarized and give rise 
to electric fields in cryogenic cavities. This effect exists 
for frequencies ranging from zero up to a few MHz. Using 
this effect all the parameters of the mirror dark matter 
plasma can in principle be measured. Besides detecting 
or constraining mirror dark matter, this experiment can 



b\ =a\ + a 2 ——-cos(S) cos(a) cos(</>2 — <j>i) , (27) 
vo 



2 Here we ignore the fact that the mirror charge density inside 

1)2 = sin(<5) sin(o!) , (28) the cavity is not constant and depends on the distance to the 

2 Vq capacitor. 



a 2 \vh\ 
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also significantly improve the upper limit on the photon 
mass. 

If the DAMA annual modulation signal is due to mirror 
dark matter then the expected signal strength, given by 
Eq.El is very small and a large cavity is therefore needed. 
The main difficulty in such a large scale experiment will 
be to keep the cavity at cryogenic temperatures. The 
required cryogenics facilities would have to be similar to 
those for the magnets of the ATLAS experiment at the 
LHC H3- 

There are however scenarios which can yield a positive 
result in less sensitive, smaller scale experiments. From 
Eq. we see that the signal strength is proportional to 
e 4 , so mirror matter scenarios with a larger value for 
|e| than the favored value of |e| ~ 5 x 1CP 9 would yield 
much stronger signals. The BBN bound on |e| is 3 x 10 -8 
|29|. but increasing |e| to this value would require lower 
abundances of elements heavier than helium by a factor 
of about 100 to avoid conflict with the DAMA results. 
This is problematic because we need to assume a high 
rate of supernova explosions to keep the mirror plasma 
hot enough to prevent it from collapsing into a disk |55j . 

It may be possible to evade these difficulties in broken 
mirror matter models Q. Here each mirror particle is 
heavier than its ordinary partner by some constant fac- 
tor. This causes the rate of cooling of the halo to be 
slower than in the unbroken case; the time needed for 
the halo to collapse is ~ 10 8 years times the mass ratio 
squared |56| . If a photon-mirror photon kinetic mixing 
also exists in this case then this type of dark matter can 
be detected in electromagnetic field penetration exper- 
iments. The signal can be much stronger than in the 
unbroken case if conflicts with dark matter direct detec- 
tion experiments can be avoided. This is possible if the 
mirror elements are not heavier than about 15 atomic 
mass units. As explained in Ref. |5jfl , the mirror neutron 
can be unstable, even if bound in nuclei. As a conse- 
quence, mirror hydrogen would then be the only stable 
mirror element. This means that if the mirror masses are 
10 to 15 times larger then one can have both a spherical 
halo and no signal in dark matter direct detection exper- 
iments. The DAMA signal would then have to be due to 
a different component of the dark matter or not due to 
dark matter at all. Because mirror electrons are heav- 
ier in this case, the BBN bound on |e| can be relaxed a 
bit |58| . Also there is no positronium- mirror positronium 
oscillation in this case. 

But whatever the theoretical motivations are for doing 
this experiment, the fact remains that the strongest null 
result to date comes from the most precise direct test 
of Coulomb's law [5^], which corresponds to an effective 
upper limit on |e| of2x 10~ 5 . This means that unknown 
physics could yield signal strengths of up to ~ 10 14 times 
the estimate given in Eq. |5J 
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APPENDIX A: SOLVING THE VLASOV 
EQUATION 

In this section we'll solve the Vlasov equation to first 
order in perturbation theory. The methods used here can 
be found in books on plasma physics, such as e.g. Refs. 
[48l l53l | . In section 1 1 i 1 1 we obtained the first order result 
03: 

^+v-Vh~—(E + vxB)- V„/o = 0, 
ot m c 

V-E = eV-E e *-- [ d 3 vf u (A1) 
£0 J 

B = ei3ext, 

where /o is the unperturbed mirror electron distribution 
given by Eq. 1131 and /1 is the first order correction to 
the mirror electron distribution. Attempting to solve 
these equations by performing a Fourier transformation 
w.r.t. the space and time coordinates leads to singulari- 
ties. This is caused by the absence of dissipative terms 
in the Vlasov equations. To deal with this problem one 
must treat the Vlasov equations as an initial value prob- 
lem. This is most conveniently done by performing a 
Fourier transformation w.r.t. space and a Laplace trans- 
formation w.r.t. time. In the following we shall use the 
notations: 



F(k,p)= dtF(k,t)exp(-pt). 



(A2) 



We now assume that at t = 0, /1 and the external fields 
are zero. Applying a Fourier-Laplace transformation to 
(|aT)| gives: 



Ph +iv-kf x (E + v x B) • V„/ = 0, 

ik-E^eik- E cxt - — [ d 3 vh, ( A3 ) 
£0 J 

B = eB ext . 

Eliminating /1 from these equations yields the mirror 
charge density in terms of the external fields: 



-_ - 1 _e ff 

P — e Pcxt + e 2X^k 2 ^ Pexk- 



(A4) 



i(k, P ) 

Here we have defined I(k,p) as: 

I(k,p) = -^= [ dx ex v(~ x ) (A -, 

\/7T ./_ 



( t — *P 4- ^h-fc 
[ X V Q \k\ ^ u |fc| 



8 



Im(p) 



3) 




Re(p) 



FIG. 3: Deforming the Bromwich contour. The contour on 
the right is the Bromwich contour C in Eq. IA9I the dots are 
the poles and the deformed contour is shown on the left. 



and p^t is defined as 



Fcxt 



e V ■ [vh x B cxt + E ext ] 



We now assume that p° xt is of the form 
p^ t (x,t) = A(x) sin(cjt). 
This means that p^xt 1S given by: 



2i 



1 



1 



p — iui p + iui 



(A6) 
(A7) 

(A8) 



Inserting this in i|A4(l and taking the inverse Laplace 
transform gives: 



P = epext + 



A(k) 



2m 2i 



dp 



exp(pt) 

2 A 2 fc 2 



J(fc,p) 



- 1 



1 



1 



p — iuj p + iuj 



r ' A!n 

where C is a Bromwich contour with Re(p) > 0. We want 
to find p at large times when the transient effects, caused 
by the switching on of the external fields, have vanished. 
The standard technique is to deform the Bromwich con- 
tour by moving it to the left of the imaginary axis but 
such that the poles are still encircled to the right, see 
Fig- El Then, for t — > oo, only the contribution from the 
poles at p = ±iu> will survive. Because the integral defin- 
ing I(k,p) in Eg. IA5I is singular for Re(p) = we must 
analytically continue it from Re(p) > to the region 
Re(p) < 0. This can be done by deforming the integra- 
tion contour in Eq. IA5I by letting it pass below the pole 
at x — — ~^7t — ^^nr . Denoting the analytical continuation 

v \k\ v \k\ o •> 

of I(k,p) by I*(k 1 p) we can then write: 



P = epext+e 



2i 



exp(io;t) exp(— iojt) 



2A 2 3 fc 2 



2A^fc 2 
I*(k,—icj) 



(A10) 



An expression for J*(fc,p) in terms of error functions can 
be found as follows. It is convenient to introduce the 
plasma dispersion function Z(Q which is defined as: 



HO 



1- [ d Z C ^~ z2 \ 

A Jc- z ~C 



(All) 



where C~ is a contour from — oo to +oo that passes be- 
low the pole at z = £. The function I*(k,p) can then be 
expressed in terms of the derivative of the plasma disper- 
sion function as 



r(k, P ) = z' 



ip 



Vh ■ k 



y v \k\ v \k\ 
Applying a partial integration to 



z'(0 



i 



dz 



exp(— z 2 ) 



'* Jc- (z-CY 
leads to the differential equation: 

Z'(C) + 2(^(0 + 2 = 0. 



(A12) 



(A13) 



(A14) 



Solving this equation and using that Z(0) = i^pH gives: 



Z(C) = V^exp(-C 2 )[*-crfi(C)] 
where we have defined: 



erfi(C) 



erf«) 



C 



exp(i ) dt. 



We thus see that Z'(() is given by: 

Z'(Q = -2 + 2^Q ex P (-C 2 ) (erfi(C) - i) • 



(A15) 



(A16) 



(A17) 



Finally, applying an inverse Fourier transformation to 
(|A10fl gives the result: 

p{x 7 t) =ep oxt (x,i) + 

- / d 3 yA(y)[G(x - y,u,t) - G(x - y,-u,t)] , 

(A18) 



where 

G{x,LU,t) 



o- / d kexpUk ■ x) X 

2i(27r) 3 ' 



2X 2 D k 2 



LP __ Vh_-k 

V^uJ-pXulkl v \k\ 



z' - 



n -1 



(A19) 



APPENDIX B: TIME DEPENDENCE OF COS(6») 

The values for 5 and t in Eq. 1251 can be obtained from 
the angle of 60° the ecliptic plane tilts relative to the 



9 



normal of the galactic equatorial plane in the direction 
of motion of the Sun around the Milky Way, the date 
of 2 July when the component of the Earth's velocity in 
the direction of the motion of the Sun around the Milky 
Way is maximal (at that time the Earth moves toward 
the galactic South Pole) an d the fact that the Earth's 
axis is tilted by 23.5° relative to the normal of the ecliptic 
plane. Expressing cos(<5) as the inner product of the unit 
vectors pointing in the directions of — Vh and the Earth's 
axis gives: 

cos(<5) = cos(30°) cos(23.5°) + sin(30°) sin(23.5°) cos( 7 ) . 

(Bl) 

Here 7 is the angle between — Vh projected onto the eclip- 
tic plane and the Earth's axis projected onto this plane. 
Approximating Earth's orbit as a circle, and using that 
on 21 June the Earth's axis points toward the Sun while 
on 2nd June —Vh projected onto the ecliptic is parallel to 
Earth's orbit, we find that 

7 » 90° + 19 dayS 360° « 108.7°. (B2) 
1 year 



Inserting this in i|Bl|) gives 6 w 43°. 

To compute to, we note that around 21 March when the 
Earth's axis is orthogonal to the unit vector s pointing 
from the Earth to the Sun, hour local sidereal time 
corresponds to the moment when s projected onto the 
Earth's equatorial plane and the position of the detector 
projected onto that plane are parallel. The angle between 
s projected onto the Earth's equatorial plane and — i>h 
projected onto this plane expressed in hours thus gives 
U). This angle follows from the inner product of s with 
— Vh by using the fact that s lies in the intersection of the 
ecliptic plane and the Earth's equatorial plane: 



-Vh ■ s = cos(60°) cos(7 - 90°) = sm(6) cos(t ) ■ (B3) 



This, combined with the fact that to must be in the in- 
terval ranging from 270° and 360°, gives to « 314° w 21 
hours. 
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